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The magneto-elastic interaction in cubic helimagnets with B20 symmetry is considered. It is shown 
that this interaction is responsible for negative contribution to the square of the spin-wave gap A 
which is alone has to disrupt assumed helical structure. It is suggested that competition between 
positive part of Af which stems from magnon-magnon interaction and its negative magneto-elastic 
part leads to the quantum phase transition observed at high pressure in MnSi and FeGe. This 
00 ' transition has to occur when A^ — 0. For MnSi from rough estimations at ambient pressure both 

' parts Aj and |Aa/_b| are comparable with the experimentally observed gap. The magneto-elastic 

\ interaction is responsible also for 2k modulation of the lattice where k is the helix wave-vector and 

OA ■ contribution to the magnetic anisotropy. 

}_( ' Experimental observation by x-ray and neutron scattering the lattice modulation allows determine 

^ l| the strength of anisotropic part of the magneto-elastic interaction responsible for above phenomena 

. and the lattice helicity. 

PACS numbers: 05.70.F h, 75.20.E n, 75.45+j 



I. INTRODUCTION 

Helical magnetic order stipulated by Dzyaloshinskii-Moriya interaction and its evolution in magnetic field attract 
■ a lot of attention in last years (see for example^ and references therein). In this respect non-centrosymmetric cubic 
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helimagnets such as MnSi, FeGe, FeCoSi with P2i3 {B20) symmetry play special role. They are a subject of 
the intensive experimental and theoretical studies during several decades. Their single-handed helical structure was 
explained by Dzyaloshinskii^. The full set of interactions responsible for observed helical structure (Bak- Jensen model) 
, was established later in^^ in agreement with existing experimental data (see for example*^ and references therein). 
(~| ' The renascence in this field began with a discovery of the quantum phase transition to a disordered (partially ordered) 
O , state in MnSi at high pressure (magnetization and resistivity measurements^ji and neutron scattering2,j2,) . Recently 
, ^ , ' similar transition was observed in FeGe^. The following properties of this state attract the main attention: i) non- 
[ Fermi-liquid conductivity, ii) spherical neutron scattering surface with the weak maxima along the (110) axea^fii, 
T-H ■ whereas at ambient pressure Bragg reflections and critical scattering were observed along (111)-^. These features and 
K*" ' the structure of the partially ordered state were discussed in several theoretical papers (se o^^i^^'^'^i^^ and references 
00 , therein). It should be noted that the spherical scattering surface with maxima along (111) was observed at ambient 
' pressure just above critical temperature Tc ~ 29K and explained using the Bak- Jensen model^^. It was demonstrated 
, also that this phase transition is of the first order oneii Recently phase separation in MnSi near the quantum phase 
• ' transition was observed by muon spin resonant method^. 
I IniSil strong anisotropy of the spin- wave spectrum at low momenta was demonstrated: excitations with momentum 
along and perpendicular to the helix wave-vector k have linear and quadratic dispersion respectively. At the same 
^ time there is a contradiction between these papers. Ir*i^ was claimed that the spin-waves are gapless Goldstone 
• • , excitations due to translation invariance along the helix axis whereas ini^ the spin- wave gap was calculated in l/S* 
approximation. This contradiction was discussed in^°. In brief its essence is following. In^** the l/S corrections to the 
spin-wave energy were not evaluated and the translation invariance was not proved. Meanwhile doing this the authors 
jj] ■ should meet a problem how to consider the Dzyaloshinskii-Moriya interaction. It contains two spin operators and if 
^ ' they belong to the single lattice point, the translation invariance holds and the gap is zero. However this interaction 
always acts between different spins, the translation changes the pair energy and the gap is not zero^i:. 

Existence of the gap A is very important for correct description of the helix behavior in magnetic field H± per- 
pendicular to the helix vector k. In the gapless case the spin-wave spectrum becomes unstable in infinitesimal H± in 
contradiction to the well-known experimental findings^^ and predictions of the phcnomcnological Landau- like theory^^. 
Ini^ was shown that the helix state remains stable if H± < and then k begins rotate toward the field. Recently 

this prediction was confirmed using small angle polarized neutron scattering and was found that A ~ 13fieV for 

The quantum phase transition to magnetically disordered state observed in MnSi and i^eCe^i^ '^^i^^ is an another 
important problem. Recently it was considered on the base of phenomenological Landau theory^!. However up 
to now we have not any attempt to understand microscopic origin of this transition. In this paper we consider the 
magneto-elastic (ME) interaction and evaluate its contribution to the square of the spin- wave gap A. We demonstrate 
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that 

A2 = A? + AL^, (1) 

where Aj was evaluated ir^ and A^^^ < appears due to the ME considered as the second order perturbation. It 
is important to note that if Aj is zero the ME has to disrupt the hehcal magnetic order. 

Rough estimations using existing experimental data at ambient pressure (see Sec.V) give |Am_e| ~ 7.6^eV and 
A/ = 4.0 28fieV. Both contributions are comparable with experimental value determined ir^. Hence at pressure 
two parts of A^ has to compete and the quantum phase transition to the partly disordered state occurs when A = 0. 
Besides we estimated the ME contribution to the magnetic anisotropy and demonstrated that it is not small in 
comparison to the experimental value. We demonstrated also that the ME leads to the lattice deformation with the 
wave-vector 2k and evaluated intensities of corresponding super-lattice reflections. Their experimental study would 
allow determine the strength of anisotropic part of the ME interaction responsible for above mentioned phenomena 
and the lattice helicity. It has to be noted that the ME interaction in the Landau theory was investigated i n^^i^^ and 
the lattice deformation was predicted at the wave- vector k. This result does not contradict to ours as it was obtained 
for magnetized systems only. 

The paper is organized as follows. In Sec. II we consider the ME in cubic helimagnets. Classical groimd state 
energy and the lattice deformation is studied in Sec. III. The spin-wave-phonon interaction and the ME contribution 
to A^ is considered in Sec. IV. Obtained results, numerical estimations and experimental consequences are discussed in 
Sec.V and main results are summarized in Sec. IV In Appendixes A and B some mathematical details are considered. 
Appendix C is devoted to consideration super-lattice reflections near forbidden < nOO > Bragg peaks with odd n-s. 



II. MAGNETO-ELASTIC INTERACTION 



In general form the magneto elastic energy is given by (see for exampl^S) 

Vme = S^S^Baf3'yf_tUjf^{R) (2) 

R 

where is the spin components at the lattice point R , U^^ = {l/2){du^ / dR^ + du^/ dR^) is the deformation tensor 
and the lattice site displacement has well known form 

where e^j are vectors of the phonon polarization , b{b'^) their absorption (excitation) operators and c-qj = ^qj~' 
Tensor B is symmetric in (a/3) and (7/i) components. In cubic crystals we have following non-zero components^ 

-E^xxxx — Byyyy — ^zzzz — ^1: -^xyxy — Byzyz — ^zxzx — -^2; (4) 

and Bxyxy = By^xy = B^yyx etc. In isotropic medium we have 

Bap^ii, = Bis{Sap,Sf3i, + Sa^Sf3^)/2, (5) 

and Bi = B,s, B2 = B,,/2. 

In cubic helimagnets the lattice spin is given byi^ 

Sr = Ae'^-'^{S^ cos a + iS'^- sin a) + A*e~*-^(5'^ cos a-iS'^- 5^ sin a) + c{S^ sin a + cos a) 

(6) 

^AS^ + A*Si +cS^, 

where < > 7^ is an average value of the lattice spin, describe its perpendicular motion, k is the helix 
wave-vector, A = (a — ib)/2, unit vectors a, b, andc form right-handed orthogonal frame, sina — ~H\\/Hc where H\\ 
is the magnetic field component along the helix vector k and He is the critical field for transition to ferromagnetic 
state. According t(M- the vector k j] c in arbitrary field. 

Using standard definition Sq — N~^^^ ^ Sr exp (— iq • R) in momentum space we obtair^ 

Sq = 5^c + 5^A + 5f A* 

S'q = sin a + 5^ cos a, — S^_^^ cos a — S^_^^ sin a + iS'q_i^, = '5'q+k "^o^ ~ "^q+k ^^^^ ~ *'^q+k' 



3 



where S^, and are functions of q, q — k and q + k respectively. 
The spin components in Eq.(7) have wcU- known form 

Si = A^i/25Jq,o - (a+a)q; = -*V^K - «iq - (a+aX/2S]; - + «iq " (a+a')q/25], (8) 

where Oq and Cq are conventional spin-wave operators. 
In momentum space Eq.(2) is given by 

Vme^N-'/^ E ^qi^^.5'^-qi-q.' (9) 

i,m=A, A* ,c 

This expression is divided on three parts: direct (cc and A A*) terms where the U tensor is k independent, first 
{cA), (cA*) and second (A A and A* A*) order umklapp terms where U operator depends on qi + q2 ± k and 
qi + q2 ± 2k respectivelj*^^ . For their consideration we use following identity 

mnBlJQ—iBA ^ mpUpQpUp + ii32[(m • Q)(n • u) + (m • u)(n • Q)], (10) 

p—x,y,z 

where Ba = Bi — 2B2 is an anisotropic part of the tensor B and u — uq . 

In the case of uniform pressure P we have C/a/3 = —{P/'iK)Sa/3 where K is the bulk modulus and Uq ^ N^^'^Sq^- 
As a result the umklapp terms are zero asc-A = A- A = and Vme —NBiS{S + 1)P/{3K). Hence the uniform 
pressure contributes to the classical part of the magneto-elastic ground state energy only. However it has to change 
basic parameters of the problem such as -81,2, sound velocities etc. 



III. GROUND STATE ENERGY AND LATTICE DEFORMATION 



In zero magnetic field sin a = and we have a planar helix. In this case from Eqs.(7-10) follows that in the classical 
part of the ME interaction the first-order umklapps are forbidden^ and we obtain 



'2iN^^^S^BAk[{g ■ u) 



c.cj, 



(11) 



where gp = ApCpBnd Up = u^2k- 



For evaluation of the lattice deformation we must consider the elastic energy. Unlike Refs j^^i^^ for simplicity we 
ignore in it the cubic symmetry as its principal symmetry breaking role has been taken into account above in Eq(ll). 
In this case the unit cell energy is given by^ 



Fir) = Q[C/.Mr)t//3a(r) + cTUl^ir)/{l - 2a)], 



(12) 



where Q — Evo/[2{l + a)], vq is the unit cell volume, and E and a are Young modulus and Poisson coefficient 
respectively. As a result we have 



i^2k = 2Q[A:2(u • u*) + (k • u)(k • u*)/(l - 2a)]. 
From Eqs. (11-13) for the ME part of the ground state energy we obtain 

Eme = i2NS*Bl/Q){{y^ ■ w*) + (w • c)(w* • c)/(l - 2a) - [(g • w) + cc.]}, 
where u_2k = —i{BAS'^N^^'^/Qk)'w. Minimum of this energy is evaluated in Appendix A and we have 



E 



ME 



U 2k 



NS^B\{l + a) 



AEvo 

-2iN^/^S^BA{l + a) 
Evok 



(Gi — G2) + 777:; 7G2 



g 



2(1 -a) 

(g* ■c)c " 

2(1 -a)_ 



NEvogj 



(Gi - G2 + G2/2), 



i2tgA/k)[g*-{g*-c)c/2] 



(13) 
(14) 

(15) 
(16) 



where gA = S'^Ba/Evq, cubic invariants Gi = 16(g • g*) and G2 = 16(g • c)(g* • c) are considered in Appendix A and 
in the right hand side of both equations we neglect a as it is usually small^''. 
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IV. MAGNON-PHONON INTERACTION 



We consider now the magnon-phonon interaction. We are interested by terms which survive at q = and contribute 
to the spin- wave gap as other terms are small corrections to the g-dependent part of the magnon dispersion considered 
in^^ andi^. To single out them we have to replace in Eq.(9) one of S'^ operators by S'(5q±k,o- As a result we obtain 
terms with phonon momenta q, q±k, andq± 2k. The former disappear at q = 0, the second are proportional to 
flq + atq and can not contribute to the gap (see below). Using identity (10) for the last AA end A* A* terms in the 
case of planar helix we have 

V2k = -(25)3/2z/cB^ ^[(aq - aiq)ffp<2k-q + («-q ~ <).9>2k+q]. (17) 

From this equation for the magnon-magnon interaction we obtain 

Fmm = 2(2A;52B^)V^^(aq-atq) [Y, 9pDpr{A,2k)g:] (a_q-a+), (18) 
where we neglect q in comparison with ±2k . The phonon Green function can be represented as 

Dpr{uJ, Q) = DtiSpr - QpQr) + DiQpQr, (19) 

where Q — Q/Q and — [M{uj'^ — s^j)?^)]^^ where ^(t)labels longitudinal ( transverse) phonon mode and s/(t) is 
a corresponding sound velocity. We neglect optical branches as their contribution is of order [sk/Oof' ^ 1 where Od 
is Debye temperature. 

In the linear spin-wave theory the Hamiltonian is given bj^ 

Hsw - Jlt^q'^q^q + ^q(«q«-q + aiq<)/2] (20) 

and the square of the spin- wave energy = — B^. As was shown ini^ Eq = Bq ^ Ak"^ /2 where A is the spin- wave 
stiffness at q ^ k and we have gapless excitations. We assume that the ME interaction is weak and gives small 
corrections to Aq and Bq. In this case from Eqs. (18-20) for the magneto-elastic contribution to the square of the 
spin-wave gap we obtain 

Ak ^BAS^ ( G1-G2 G2 \ AeEv^ gl , 
2SM \ s? + s2 ) - 45. 



= ' ^^^^ + ^ - - (gi + g,) (21) 



where qa = S^Ba/Evq and neglecting a we have Msf = Evq and St = s;/\/2- This expression is negative as it should 
be in the second order perturbation theory and A^^^ = in < 100 > direction only (see Appendix A). Consideration 
of q± k terms lead to expression similar to Eq.(17) with replacing a — a'^ a + a'^ which do not contribute to the 
gap. 

The spin- wave interaction considered in the 1/5 approximation leads to positive contribution to which is given 

by 

and I?q is a form-factor of the Dzyaloshinskii interactioni^j^^. The helical structure can be stable if 

A^ = A? + A^,^ > (23) 

and if A/ = it can survive at k ||< 100 > only where Amb — 0. Meanwhile it is well known that in MnSi and 
FeGe at low T the helix axis k ||< 111 > and A/ > \Ame\- 



V. EXPERIMENTAL CONSEQUENCES AND DISCUSSION 



For discussion experimental consequences of the magneto-elastic interaction we have to know the Young modules 
E and anisotropic part of the magneto elastic interaction S^Ba- For MnSi according to^;^ the bulk modulus 
K = 1.37 X lO^bar and neglecting the Poisson coefficient a we obtain E = 3K = 4, 11 x 10^6ar and Evo = 2A0eV 
{vo = 95 X 10-24cm3). 
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Unfortunately the value of BaS"^ is unknown. As we will see below it may be determined by x-ray and neutron 
scattering. Isotropic part of the ME interaction was studied by indirect method in22. Its contribution to the lattice 
constant Aa/a ~ —1.1 x 10^"' at T = QK was determined and the sum of the isotropic part of the ME and elastic 
energies can be represented as 

B,,S^{Aa/a) + Kva{3Aa/a)^ /2. (24) 

This expression is minimal at gis = S'^Bis/ Evq = —3.3 x 10^^ and from Eqs.(15) and (21) we obtain 

Eme = ~6.5^ieV{57mT}{gA/g^sfiGl - G2/2), (25) 
AifB - ~{l7fieV{0.15T}f{gA/g;sfiGi+G2), (26) 

where we used S = 1.6, Ak^ = He = 0.6T and He is a critical field for transition to ferromagnetic stateii^i. 

We compare now above results with existing experimental data^'^. As was shown iniS, the ground state energy of 
the helical structure in magnetic field is given by 

SH^l SHlA? 

Eg = Ea + Eme - ^ - ^Hc{A^ ^ Hl/2y ^^^^ 

where is a field component along (perpendicular) to the helix wave- vector k, Ea = {S'^Fok'^ — 3S'^K)L/4:, Fq 

and K are constants of the anisotropic exchange and cubic anisotropy respectivelj*^^. The cubic invariant L = (g • c) 
is considered in Appendix A. 

Evolution of the helical structure in magnetic field was studied by small angle polarized neutron scattering in MnSi 
near Tc?^ at low T^^ and in compound FeCoSi^'^-^^ . Two new characteristic fields were determined. In zero field 
the sample is in multidomain state with k along all < 111 > directions. Then for the field along one of < 111 > axes 
at Hci the single domain state appears. With further field increasing the Bragg intensity demonstrated a cusp at 
Hin. In22 it was interpreted as instability of the k direction connected with the second term in Eq.(27). Indeed if H 
is slightly below A\/2 this term predominates and the vector k has to rotate perpendicular to the field but blocked 
by the anisotropy. Just below Tc where the anisotropy is weak this rotation was observed in^^. In MnSi we have 
Hci = 80mT, H,nl60mT and A ~ llOmT = 13/ieV21. 

It is obvious that the single domain state can be reahzed if SH^j/2Hc — 9mT is of order of Ea+Eme- For < 111 > 
directions we have Eme — "^^{gA/ gisY'mT [See Eq.(A5)]. So this condition is fulfilled roughly. It is impossible to do 
more detailed analysis as we do not know Ea and gA- 

The invariant L has two extrema L = 2/3 and L = at k along < 111 > and < 100 > directions respectively and 
a saddle points at < 110 > directions. Hence if one neglects the ME interaction the configuration with k ||< 110 > 
is forbidden**. The same holds for maxima of the critical fluctuations above Tc — 29 K-^. Meanwhile in MnSi at 
high pressure above the quantum critical point pc ~ lA.Gkbar maxima of the neutron scattering at < 110 > directions 
were observed*^. In Appendix B we show that the ME interaction can not resolve this problem i.e. that < 111 > and 
< 100 > remain only possible k directions in zero magnetic field. 

Let us estimate now two contribution to the spin- wave gap given by Eqs.(22) and (26). We do not know real form of 
theratior = Dq/Z^o inEq.(22). For r = l,Ak'^ = 0.6T and S* = 1.6 we get maximal value A^^^^ = (0.24r)2. Minimal 
value of Aj may be estimated assuming that in Eq.(22) qmax = 0.024nTO is a border of the Stoner continuum^. In 
this case we have Aj^j„ = (0.035T)^. In Eq.(26) for < 111 > direction expression in the bracket is equal to 4/9 [see 
Eq.(A5)] and we get Ame — 0-0Q7{gA/gis)^T . So we see that at ambient pressure both contributions are comparable 
with the observed A ~ O.llT. Hence we can do plausible assumption that the quantum phase transition observed at 
lA.Gkbar^^i^'^ is a result of vanishing A^. At higher pressure the helical structure becomes unstable. 

For more precise estimations experimental measurement of magneto-elastic anisotropy S^{Bi — 2B2) would be 
important. We demonstrate now that it can be directly extracted from intensities of satellite peaks near nuclear 
Bragg reflections. Indeed using Eq.(16) we obtain'^^ 

«±(K) - (2.g^/fc)2|K± • [g- (g • c)c/2]|2|F(K±)|2, (28) 

where K± = K±2k, F(Q) is the nuclear structure factor, K is a reciprocal lattice point and gA = S^{Bi — 2B2)/ Evq. 
Relative satellite intensities are given by 

5J±(K)//(K) ^ {2gA/k)'\K • [g - (g • c)c/2]|2. (29) 

In zero magnetic fleld vectors k are along all < 111 > directions. If K || k we have 5I± = 0. If however k || (1, 1, 1) 
but K = (27rn/a)(l, 1,-1) we obtain 



, ^ , = 2.2 X 10-^ ^ (30) 
\9V3kaJ V 5*. 
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where we have used ka = 0.17~. For K || (1,±1,0) this expression has to be multiphed on 1/4. 

Similar results can be obtained for other K directions with one exclusion. We are interested by crystals with 
P2i3 {B20) symmetry where < n, 0,0 > Bragg reflections are forbidden if n is odd and observation of very weak 
super-lattice reflections would be more easier^*'. Hence this case has to be considered separately. We restrict ourselves 
by the case k || (1, 1, 1) only. First of all for even n we have Eq.(30) with replacement 2.2 2.2/4 ~ 0.55. For the 
odd n the Bragg intensities of the satellites are given by 

/±(K) = 0.55 X 10-5 ('2Ea) |Fm„(Q) + Fs^m^ (31) 

where form-factors -Fj(Q) are given in Appendix C. They are not zero due 2k modulation. Unfortunately Eq.(31) 
has an additional small factor (2fca)^ << 1 in comparison with the even-n case. Observation of these odd reflections 
provides a possibility to determine the lattice chirality [See Eq.(C6)] and its connections with the spin chirality 
studied by polarized neutron9^^i2ii2^i2&. It has to be noted that the lattice chirality in some cases was determined 
by anomalous x-ray scattering^, and electron diffraction^^ . 

There are eight domains in virgin sample corresponding k along < 111 > directions. In magnetic field Hex — 
0.08T << He the single domain state is realized with k along the field and satellite intensity increases. However 
further increasing of the field suppresses the helical structure and it disappears at He — 0.6T— . In intermediate 
region at Hci < H < He the lattice modulation with the wave-vector k has to appears also^^. Along with discussed 
2k lattice modulation at low field H < A\/2 the second order helix harmonic appearsi^. I was observed in^'^''*^. 



VI. CONCLUSIONS 



We considered the magneto-elastic interaction in cubic helimagnets with _B20 structure and demonstrated that 
it deformed the lattice and gave a negative contribution to the square of the spin- wave gap . Hence the helical 
structure is stabilized due to positive contribution to which stems from the magnon-magnon interaction^^. It 
was suggested that the quantum phase transition observed at pressure in MnSi and FeGe is a result o competition 
between these two parts of the gap and takes place when A^ = 0. This suggestion is supported by rough estimations 
at ambient pressure of both contributions to A^for MnSi which have the same order and close to experimentally 
observed gap. It was discussed also how to measure directly anisotropic part of the ME interaction responsible for 
considered phenomena using x-ray and neutron scattering. 
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APPENDIX A 



In this Appendix we calculate deformation of the lattice by the ME interactions, cubic invariants G1.2 in Eqs.(15), 
(21), (25-26) and analyse their properties. 
We begin with the classical energy (14). It is minimal if 

Wp + {^- c)cp/(l - 2ct) = 5;, (Al) 

where p = x,y,z and Qp — A^Cp and we have 

Wp ^ 9*- (c- g*)cp/[2(l - cr)], 
EcL = -{2NS'ByQ)[ig ■ g*) - (g . c)(c • g*)/(l - a)] (A2) 

where 16(g • g*) = Gi and 16(g • c)(c • g*) = G2. 
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In cubic xyz frame we can write 

a = (cos I? cos if, COS'S sin ip, — sin -d); b — (sin yj, — cos </?, 0) ; c — (sin z9 cos (/s, sin i? sin cos , ( A-3) 

and for three principal k-directions < 111 >, < 110 >,and < 100 > we have: {a — {l/\/6, 1/V6, ~\/2/3); b = 

(1,-1,0)/V2; c = (1,1,1)/V3}, {a = (0,0,-1); S = (1,-1,0)/V2; c = (1,1,0)/V2}, {a = (0,0,-1); 6 = 
(0,-1,0); c(l,0,0)} respectively. 
In this representation for G-functions we obtain 

Gi = sm^i?[(cos^ cos^ + sin^ (^)^ cos^ + (cos^ ■& sin^ <y9 + cos^ ip)'^ sin^ <y9 + sin^ cos^ i9] ^^^^ 
6*2 = sin^ ')?{[cos^ 'i?(l + sin^ ip + cos"* if) + 2 sin^ ip cos^ ip]"^ + 4(sin'^ — cos^ ip)"^ cos^ sin^ ip cos^ (/j} 

From these equations follows that functions Gi — G2 and G2 have extrema at<lll>, <110> and < 100 > near 
which we have 

Gi - G2 = 4/9 - 20(5i?V9 - 40V/27; G2 = S6^^/9 + 16^/27, < 111 >; 

Gi - G2 = m^^/i + Sip^; G2 = 1/4 - 2{6'd^ + V), < HO >; (A5) 
Gi - G2 = S-d^ + Sip'^; G2 = iiS-d"^ + Sip-^ + Sd^Sip^), < 100 >, 

where S'd and 6ip are distances from corresponding extremal points. Hence in considered directions both functions 
Gi — G2 and G2 have extrema and one can show that they have not other extrema. 

Contribution of the anisotropic exchange and cubic anisotropy to the classical energy is proportional to cubic 
invariant L given by^^ 

L = 4^ l^pPcp = sin^ i?[(cos^ ?9cos^ ip + sin^ (p) cos^ p + (cos^ i9sin^ p + cos^ p) sin^ (p + cos^ ■&]. (A6) 
As above for three principal directions we have 

L = 2/3 - 4(5z9V3 - S6ip^/9, < 111 >; 1/2 + S^"^ - 25lp^ , < 110 >; 2{Sd'^ + 6ip^), < 100 >, (A7) 
and L has a saddle point at < 110 >. 



APPENDIX B 

We demonstrate now that in presence of the ME contribution to the ground state energy the < 111 > and < 100 > 
remain only possible stable directions for the vector k. 
From Eq.(27) at H = follows 

£;g = $L - *(Gi - G2 + G2/2) = */(y), (Bl) 

where ^ff > and y = ^/^f. 

We have to study behavior of f{y) for three principal directions. For k |1< 111 > we obtain 

fiy) = (2/3)(y - 2/3) + (4/3)(-j/ + 4/3)iSd' + 2Sp^). (B2) 

and Eg is stable if <i> < 4^/3. In the < 110 > case we have 

f{y) = (l/2)(2; - 1/4) + {y - 9/i)S{)' - 2y V, (B3) 

In this configuration there is a saddle point as coefficients at deviations (5t?^ and Spi'^ can not be positive simultaneously. 
Finally if k ||< 100 > we have 

f{y) = 2{y~l/2){Si)^ + S^'). (B4) 

This configuration is stable if $ > '5/2. However, comparing Eqs.(B2) and (B4) we see that configuration < 111 > 
in the region 5'/2 < $ < ^'/3 has lower energy and < 100 > configuration is metastable. Hence we see that the 
magneto-elastic energy can not be responsible for the stability of < 110 > configuration. 
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APPENDIX C 

There are two different ions in compounds with P2i3 symmetry {Mn and S'l'^'S.; Fe and etc.) labeled 

below as 1 and 2 respectively. Each of them occupy in cubic unit cell four positions: pi = {x,x,x), p2 — (1/2 + 
x, 1/2 — x, 1 — x), p3 = (1 — a;, 1/2 + a;, 1/2 — x) and p4 = (1/2 — 1 — x, 1/2 + x) (Right-handed structure) or 
pi = {x,x,x) p2 = (1/2 - X, 1/2 + x, 1 - x), (03 = (1/2 + 1 - 1/2 - x), and p4 = (1 - 1/2 - x, 1/2 + x) (Left- 
handed structure)^.. For MnSi we have Xi = 0.138 {Mn) and 2:2 = 0.846 {Si). It is interesting to note that these 
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From Eqs.(C4,5) follow the intensity ratio 

/j?(Q3) ^ /l(Q4) 
/fl(Q4) /l(Q3)' 
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